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r^ . Abstract 

A geometric formulation of the linear beam dynamics in accelerator physics 
is presented. In particular, it is proved that the linear transverse and longi- 
tudinal dynamics can be interpret geometrically as an approximation to the 
Jacobi equation of an affine averaged Lorentz connection. We introduce a spe- 
<~i \ cific notion reference trajectory as integral curves of the main velocity vector 

field. A perturbation caused by the statistical nature of the bunch of particles 
is considered. 



1 Introduction 
> 

jy-v . The notion of reference trajectory is fundamental in particle accelerator beam dy- 

f^ ' namics jlO^ Vol. I]. The reference trajectory is a curve A : I — )• R^ associated to the 

VO ■ motion of a bunch of particles in a particle accelerator. The reference trajectory is 

^r^ \ ideally the trajectory that one point charged particle will follow under the specifi- 

^) • cations of the initial conditions and external fields. The beam dynamics is adapted 

CN I to the reference trajectory in such a way that the dynamics of particles composing 

each bunch are controlled by electric and magnetic fields, theoretically along trajec- 
tories near to the reference trajectory. Indeed, beam parameters like emittance or 
dispersion and some observable quantities are referred to the motion relative to the 
reference trajectory. 

However, there is no theoretical guarantee that one can associate to the reference 
trajectory a real trajectory of a charged point particle or the real collective motion 
of the bunch. This is because the observation of the motion of the bunch is indirect. 
Therefore, there is the assumption that the reference trajectory is related with a 
physical trajectory of the bunch of particles, although one does not have theoretical 
proof of that fact. 

In this paper we introduce a theory where the reference trajectory is linked with 
an observable quantity that satisfies a geodesic differential equation. In particular, 
given a bunch of particles described by the kinetic function f{x,y) [H [21 [3], we 
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consider the averaged Lorentz dynamics and the corresponding geodesic flow [5]. 
It determines trajectories associated with the averaged velocity field of the kinetic 
model, that we promote to be the expected velocity field. 

The theory contains two different kind of elements. The first one is the one-particle 
distribution, that one takes as give from kinetic model. Although we will assume that 
f{x,y) is a solution of the Maxwell- Vlasov system or easier, the Vlasov's equation, 
the present paper only assumes that the distribution exists and is (weakly) smooth. 
We assume that such distribution function is given. The second kind of elements 
are from geometric nature, and related with the geometric theory of differential 
equations. In this paper we concentrate in the second aspect. We present two 
results. The first is that under some conditions, the standard linear accelerator 
particle beam dynamics corresponds to the Jacobi equation of the averaged Lorentz 
connection. Secondly, we present an observable which depends on the collective 
nature of the bunch of particles. 

This work is organized as follows. In section 2 we introduce the Lorentz connec- 
tion and the averaged Lorentz connection. In section 3 we discuss briefly the Jacobi 
equation of the averaged Lorentz connection. In section 4 we show that the linear 
transversal dynamics in accelerator physics corresponds with a convenient approxi- 
mation for the Jacobi equation of the averaged Lorentz connection. In section 5 we 
flnd an observable which is related with the composed nature of the bunch of parti- 
cles in a beam. Section 6 is devoted to the longitudinal dynamics from the Jacobi 
equation point of view. Finally, a short discussion and perspectives for future work 
are briefly commented. 

2 The Lorentz connection and the averaged Lorentz 
connection 

After introducing some notation, we define the Lorentz force connection and the 
averaged Lorentz connection and equations. We describe how the Lorentz force 
equation and the averaged Lorentz force equation are related. 

2.1 Linear connections on pull-back bundles 

Let (M, 77) be the Minkowski space-time, TM the tangent bundle, N = TM \ {0} 
and X) the unit tangent hyperboloid. The pull-back bundle 7r*TM — > N of the 
tangent bundle TM is the maximal sub-manifold of the cartesian product N x TM 
such that the following equivalence relation holds: for every u € Xl and (n, ^) G 
vrj" (u), (n, ^) G 7r*TM iff vr o tt2{u, ^) = '/r(ti), where the projection on the first and 
second factors are denoted by tti : 7r*TM — > N, (u,^) — > u, ^2 : 7r*TM — > 
TM, (n,O^C. 

A semi-spray vector field on TN, written in local coordinates on N (U, (x, y)) by 

^=G\x,y), y' = ^ (2.1) 

has associated a Hessian of the form 

r -.(x y) - i ^"^'^'^y^ (2 2) 



The proof of the fohowing result can be found in [5] , 

Proposition 2.1 Given a semi-spray G*(x, y) there is defined a connection on 7r*TM 
determined by the relations 

V^7r*Z := ^r{x, yY jk Z'' 7r*ei, Vv7r*Z := 0, V eV, Z e TTM, (2.3) 

where {7r*ej, i = 0, 1, 2, 3 } is a local frame for sections r(7r*TM). 

2.2 The Lorentz connection 

The Lorentz force equation is the second order differential equation 



d^x* „^,- dx^ dx^ ,■,• ^ dx^ I ,dx dx , , , ,, 

for t G I = [a, b], T* ^^(x) are the coefficients of the Levi-Civita connection ^V of r/ 
and F is the Faraday tensor. 

If we apply the formula ()2.2p to the semi-spray associated with the Lorentz force, 
we obtain the connection coefficients ^ 

^r ,, := ^r ,k + (F^ jly'^vmk + F' kly'^vmj) + F^ml {y'^vjk - VjsVki v'^y'v' ) • 

(2.5) 

and one can define the corresponding connection by proposition 12.11 It turns out 
that (see 0) the following proposition holds: 

Proposition 2.2 Let M, S, 7r*TM and ^V be as before. Then the auto-parallel 
curves of ^V are the solutions of the Lorentz force equation. 

2.3 The averaged Lorentz connection 

Given the Lorentz connection -^V on the bundle 7r*TM — y N we associate an 
affine connection averaging on the unit tangent hyperboloid X). In order to calculate 
average of observables, we use the measure 

1 



dvol{x,y) = ^/dO^ — dy^ A-'-dy"" \ y° = y°(x°, ...,x^y^ ...,2/^). 

The function f{x,y) is a (weak) solution of a one particle kinetic model. It can 
be associated with solutions of the Maxwell- Vlasov system or alternative dynamical 
models \XiM- Using such measure, one can prove that [5j: 

Proposition 2.3 The averaged connection of the Lorentz connection V is an affine, 
symmetric connection (^V) on M with connection coefficients 

(^r ,fe) := T ,fc + (F%-(iy'"r/,fc)r?„fc + r ,(ly'")ry„^.) + F„ i((y-) r/,fc - VjsVkliy'^y'y')). 

(2.6) 



The moments are defined by the integrals 

volCS^) := f{x,y)dvol{x,y), (?/*):= — tt^tt / y'f{x,y)dvol{x,y), 

(y'"yV):=^^ / y'^y'y'f{x,y)dvoi{x,y). 

2.4 Comparison between the geodesies of ^V and ( ^V) 

We define the energy function E of tlie one particle distribution function f{x,y) to 
be the real function 

E-.M^-R, x^ E{x) := inf{y\ y G supp{fx)}, (2.7) 

where y^ is the 0-component of the tangent velocity vector, measured in the labora- 
tory coordinate frame. The function a{x) is the diameter of the support of fx] the 
a is defined as a = supp{a{x)^ x € M}. The parameter s is the coordinate time 
measured in the laboratory frame. 

The relation between the geodesies of ^V and the geodesies of (^V) by the fol- 
lowing result O [7j 

Theorem 2.4 Let (M, 77, [A\) he a semi-Randers space. Let us assume that 

1. The auto-parallel curves of unit velocity of the connections ^V and (^V) are 
defined ont^ I. 

2. The ultra-relativistic limit holds: the energy function of the beam E{x{s)) is 
much larger than the rest mass of the particles, E{x{s))))l for all t € I. 

3. The distribution function is narrow f{x,y) = in the sense that a « 1 for 
all X € M. 

4- The support of the distribution function f{x,y) is invariant under the flow of 
the Lorentz force equation. 

5. The change in the energy function is adiabatic in the sense that -^logE{x{s)) « 
1. 

Then for the same arbitrary initial condition (x(0),2;(0)), the solutions of the equa- 
tions 

^ViX = 0, (^V)^i = 

differ in such a way that 

\\x{s)- x{s)\\ < 2(C7(x)||F||(x) + E-\s) Cl{x){l + B2{x)a))a'^ E~^{s) s'^ , (2.8) 

where the functions C(x), C2(x) and B2{x) are bounded by constants of order 1 and 
the distances are measured in the laboratory frame. 

One of the consequences of this result is that in the ultra-relativistic limit, and for 
narrow distributions, the solutions of the Lorentz force equation can be approximate 
by geodesies of the averaged Lorentz connection (^V), 

(^Vi)i = 0. (2.9) 



The integral curves of the mean velocity field V{x) ~ {y)ix) are approximate solu- 
tions of equation 12.91 (see for instance [6j ) , 

\\{'^V)yV{x)\\<C{K).a'+0{a^). (2.10) 

where C(K ) is a constant. Therefore, the geodesies of V can be approximated by 
the integral curves of the mean vector field {y){x) as equation ()2.10p showq^- 

3 The Jacobi equation of the averaged Lorentz dynam- 
ics 

We consider the Jacobi equation of a general affine, torsion free connection V. The 
Jacobi equation is a linear differential equations describing the relative behavior of 
nearby geodesies. If the central geodesic is X{t) and x{t) = ^{t)+X{t) neighborhood 
geodesic is x{t) = ^(t) + X{t), then the Jacobi equation for V is the linear, ordinary 
differential equation 

V^V^{O+R{C,X)-X = 0, (3.1) 

where R{^,X) is the curvature endomorphism of the plane generated by {■?,-'^} [S]. 
The Jacobi equation can be re-written in local coordinates as 

The averaged Lorentz connection ( ^V) is an affine connection on the tangent bundle 
TM. Therefore, we can apply the standard Jacobi equation to the averaged Lorentz 
connection, 






„,r,,,,a,..ry(^^..^f). 



Prom the form of the system of differential equations ()3.3p it is clear the following, 
1. There is a term representing the inertial force: 

^,:H'r-,.,'.,^r,)(^^.2!ff)^ (3.4, 

Ax is universal, in the sense that it is independent of the particle mass. 



^In equilibrium, one makes the assumption that (y) is equal to the expected velocity field. 



2. Ai can depend on the electromagnetic field, since ^^ can depend implicitly 
on the electromagnetic field when defining the reference trajectory. The typical 
example is the reference orbit in a betatron accelerator |10j Vol. I, chapter 3]. 

4 Transversal beam dynamics from the Jacobi equation 
of the averaged Lorentz connection 

Let us consider the laboratory coordinate frame. Then a convenient global coordinate 
system {t,x^,x'^,x^) is defined as follows. The parameter t is the proper time of 
the reference trajectory, the coordinate x^ is the Euclidean length measured in the 
laboratory frame, and rc^, x^ are the orthogonal coordinates of a point in directions 
orthogonal to j^. Given the reference geodesic X : I — > M is very useful to 
consider the relative coordinates {x^ — X^,x^ — X^) respect the reference trajectory 
(transverse degrees of freedom). In this way, {C^,x'^,^^) coincide with the variables 
used in beam dynamics flO', Chapter 4]. 

4.1 Transverse dynamics and the Jacobi equation of the averaged 
Lorentz connection 

For the analysis of the transverse dynamics we assume that the difference -^ ^ 

is constant and that the external electromagnetic fields are static magnetic fields 
in Minkowski space. We will only consider the lower order terms in the degree 
a + 5 + c of the monomials C°'{-^)^^'^ appearing in the corresponding expressions, 
with^ = x{t)-X{t), e:={y)-^. 

First we linearize the equations with respect to the degree defined by the vector 
fields along the central geodesic ^, its derivatives and e. 



+ Cr-,. + 5'a,'r-,,)(i?|i^ + 2i!^f), .,,,.,,„ = 0,1,2, a 

Since e is small, we can replace {y) — > -^, obtaining the expression 
d'^C , ^dC^ dX'' fl^^i dX"^ , ^i dX"" 



dt^ dt dt \2^ ^ dt """ " dt '^' 

.ir,fl,^i dX"^ . dX"^ A dX^ dX'' 



2' ' dt ""^ ' ^ dt ""-"J dt dt 
The condition of transversal dynamics is -^ —^ is second order in ^ and its deriva- 



tives. Suppressing this term we obtain, 



dt2 ' - dt dt \2 ' dt "^V ' ^"'V2 ' dt "^V dt dt 
Therefore, at first order the differential equations are 

._d'e,dc\i ,dx^, dxux^ , dxux^ dxue 

In the transverse dynamics, if there is not dispersion one has that -;|-F* j = 0, if 
there is not dispersion. Hence, the differential equations in this regime are 

^_d\^ di^ dx^dx\. dx^dx\ , dx^de .... 

The last term in ()4.ip is an inertial term. We assume that ^ is small compared with 
the curvature radius p of the central geodesic. Then the inertial terms are 

The last term in the fourth equation corresponds to the relative centripetal force. 

4.2 Examples of transverse linear dynamics 

We study some examples of transverse linear dynamics using the linearized version 
of the averaged Jacobi equation. 

1. Motion in a normal magnetic dipole 

The reference frame is the laboratory frame. In this case the electromagnetic 
field is given by the expression 

/ 

&o 

-6o 

V 

where 6o is the dipole strength. Since the magnetic field is constant, ^ diF^ j = 
0. Therefore, the equations of motion for the transverse degrees of freedom 
(e\aare 

dt'' ^^ dt> V' dt'' 




One can change the parameter of the curves by t — )• x^ = /. Then one has that 

— — + — = — ^ = 

These are the standard equations for a normal dipole. 

2. Motion in a normal quadrupole field combined with a dipole 

In this case the electromagnetic field has the form 



F(x) 



The Jacobi equation reduces to 

Let us consider the contributions from -^£}diF^ j and -^^^diF^ j. Using Eu- 
ler's theorem on homogenous functions one gets the relations: 

Then the differential equations are 

These are the equations of the linear transverse dynamics in quadrupoles com- 
bined with magnetic dipole field. The solutions are parameterized by the 
proper time t. If we use the Euclidean length Z, the equations are 

^-^''•.^"O, ^ + «'''. = »^ (4.2) 

These equations coincide with the standard equations in beam dynamics |10^ 
Vol.1, chapter 4]. 

3. Motion in a normal dipole combined with a 45 degrees quadrupole 

In this case the electromagnetic field is 



6o + hii^ 

-6o - he bie 

-hi^ 

Following the same procedure as before, we get the Jacobi equations 

^ + «''.4 = ». -ir-e»« = o^ (4.3) 



The above examples show how the hnear transverse dynamics can be obtained from 
the Jacobi equation of the averaged Lorentz connection. 

5 A perturbative collective dynamical effect 

We will consider the perturbation in the equation of a point charged particle caused 
by the fact that the bunch of particles is composed by a collection of many identical 
particles. 

5.1 Calculation of the dispersion function in beam dynamics 

We follow the formalism developed in flOi Chapter 4] for the treatment of linear per- 
turbations and dispersion. However, we maintain the proper time t as the parameter 
of the curves. 

As we showed in the preceding section, the transverse dynamics is determined by 
equations of the form 

^+K{t)u = (5.1) 

The general solution of ()5.ip is 

, , ^ , , ^ , , duo du , , dC , , dS , . dur, 

uit) = Cit)uo + Sit)^, -it) = -(t)no + ^(t)^, 

with initial conditions 

C(0) = 1, ^(0) = 0; S{0) = 0, ^(0) = 1 

for arbitrary initial values uq and -^. The functions C{t) and S{t) satisfy 

-^(t) + K{t)Sit) = 0, ^(t) + K{t)S{t) = 0. 

The perturbed equation has the form: 

d'^u 

— (t)+K{t)u{t)p{t). (5.2) 

A particular solution for (j5.2|) is 

P{t) = [ p{i)G{t, i)di, (5.3) 

Jo 

where G{t, t) is the Green function associated to the differential equation (j5.2p . One 
can prove that in the absence of dissipative forces (that is, which do not depend on 
the velocity of the particle), the Green function of the differential equation is 

G{t, i) = S{t)C{i) - C{t)S{i). (5.4) 



Therefore, the general solution for the equation ()5.2p is 

u{t) = au C{t) + bu S{t) + P{t). (5.5) 

This solution breaks down if there are strong loses of energy by synchrotron radiation 
or other dissipative effects. 

If the components of a bunch do not have the same energy, one obtains for the 
transverse degrees of freedom the following differential equation [10' pg 109], 



+ K{t)u = -{t)5u, A = ^, 5p = ViSpiY + (<5p2)2 + (<5p3)2. 
dt'^ po po 

We assign to 6p the maximal value of {||.^(X(t))||}, using the Euclidean metric 
defined on the laboratory coordinate system. Then the general solution for u is 
linear in the perturbation, 

u{t) = au x{t) + bu S{t) + A D{t) := a„ x{t) + 6„ S{t) + Offu{t). 

where a^ and bu are constants that depend on the initial values and D{t) is the 
dispersion function as in |101 Chapter 4]. 

5.2 An observable associated with the composed nature of the 
bunch of particles 

We have shown in the previous section that at first order and when we take the 
approximation (y) — > -^, the differential equation for the transverse motion is the 
Jacobi equation of the averaged connection. Therefore we can consider the terms on 
e^ in the averaged Jacobi equation as a perturbation and apply the method of the 
Green function. 

The off-set function is defined as 

Offl'^t) = ni'3(t) - ai,3 C'^it) - 6i,3 5i'3(t), 

the super- index refers to the transverse directions (^^ or ^^ in the laboratory 
frame defined previously. Using the corresponding Green function we obtain 

Officii) = [ p'^^it)Guit,i)di. (5.6) 

Jo 

The perturbation p^'^{t) is in this case defined by all the terms of the averaged 
Jacobi equation which are not contained in the linearized equation respect to the 
degree {a + b + c). Therefore let us re- write the Jacobi equation of the averaged 
connection. Using e^ = (y^) -if- we get 



Offl'\t)= [ dt2^(t).^(i(Fi'^(t)6,(t)+Fi'^(t)6,(t)) 
Jo "t dt ^^ 

+ ^^^(Fi'3^(t) ( (y"^)(t)ry,, - (y™y'^y'=)(t)r/,aryzfc) + 
^tdi (Fi'3^(t) ((y-)(t)r?,, - (y'"yV)(t)r/,ar/;fc))y (5.7) 
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In the integrand of equation (j5.7p we can make the substitution 

We can also consider the derivatives in the longitudinal and temporal directions. 
With a convenient choice of the coefficients 09,2 and 60,2- Then we can write 

^(t)^(a,-(t) + b,-ii)+—^), ,=0,1,2,3. 

In this expression repeated indices are not summed. For the transverse degrees of 
freedom, the unperturbed solutions are the same as before. 

For the longitudinal j = 2 and temporal j = degrees of freedom, one gets the 
following relations by comparison with the Jacobi equation, 

d^2 , , (iC72 d52 Of flit). 

Let us consider the regime where Off? = Off? = 0, Vn. Then 

dC\~, , dS\~, ^..,v~^A.~^ dX^f^/^\ 

j=0 



Offl'Ht) = 1^ di ( J2 2(«. ^' (*~) + b, §' (t) + Offiii)') (i) . ^ (i ( Fi'3 , (t) e,{i) 



;il,3 /r\ , /r\\ I dX dX /„i 3 (2\ / I„,m\i2\^ /„,m„,a„,k\ 



+ (a, C'(t) + bi 5'(t) + Offi{i))di (Fi'3„ ( (y™)(t)r?,fe - (y'"y'^y')(t)r?,ar/ife))). 

This is an integro-differential equation for Offc' ■ In the Born approximation one 
puts OffUi) = in the integrand, 

Offl'Ht) = 1^ £{2^ ■ ^{l{F''',ii)e,ii) +F'\(t)e,(t)) 

+edi (F^'3^(t) ( {vmtHk - {y"'y''y'mr^,avik))). (5.8) 

If the perturbation is constant along the reference trajectory, the term containing 
derivatives are neglected and one obtains 

0/4''(i) = J^ di{(2J^^^ ■ (Fi-%(t)e,(t) +Fi-^(t)6,(t))) + 
dX^ dX^ / \\ 

+ -^^^[F^'^mii) ( {yn{i)vjk - {y"'y''y'){i)rijamk))). (5.9) 
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This expression depends on the particular solution -^{t). One way to eliminate such 
dependence is to consider the average of the above expression with y = -^ for all 
possible geodesies (note that then ^^ = e^), 

+ ^^(f'''-W ( {ymiHk - {y'^y^y'mv.avik)) 
However, the average of the first integral is zero. Therefore. 

{Offl''){t) = j^ di[-^^^[F'''Ui){{yn{i)rijk-{y"'y''y'){iHamk)). 

(5.10) 

In the case of a delta function distribution one has that {Of ft' ){t) = 0. Therefore, 

the averaged off-set function {Of fc' ){t) is a measurement of which good is the 
approximation of the standard beam dynamics as description of the dynamics of a 
system of particles. 

The averaged off-set {Of fp' ){t) is determined by: 

1. The reference trajectory X{t), which is a geodesic of the averaged connection. 
This is theoretically fixed and corresponds to an observable quantity. 

2. The external electromagnetic field F^'^m(X(t)) along the reference trajectory. 

3. The first and third momentum moments of the distribution function /(x, y) 
along the reference trajectory. 

This properties make the function {Of fc' ){t) to be an observable associated with 
the beam of particles. 

There is a dependence on how the averages are calculated. This depends on 
the specific kinetic model that we are considering. One this is fixed, the function 
{Offc' ){t) is useful for two things, depending on the particular situation: if the 
kinetic model is well known, it can be used to test the dynamics of the accelerator; 
if the dynamics of the beam is well known, it can used to test the kinetic model 
used. It is also interesting that indeed it depends on only three moments of the 
distribution function. 



6 Longitudinal beam dynamics and corrections from 
the Jacobi equation of the averaged Lorentz connec- 
tion 

Longitudinal dynamics is the theory in beam dynamics that deal with the accelera- 
tion mechanisms [lUt Vol. I, chapter 8]. We show that in both, the Jacobi equation 
of the averaged Lorentz dynamics provides equations for longitudinal dynamics, even 
if in an indirect way. Two situations are considered: a constant electric field and an 
alternating electric field. 

12 



Let us consider an inertial coordinate system defined by the vector field Z = -^ 
corresponding to tlie laboratory frame. The interaction of an ultra-relativistic bunch 
of particles with an external longitudinal electric field E = {0,E2{x),0) and zero 

magnetic field. For narrow distributions it holds that -^-^ = O^. This relation 
can be seen as follows. For the linear dynamics S, = (^, 0, —S,, 0) in the laboratory 

frame. Also note that in the ultra-relativistic limit ^^ = (1 + E,0,E,0), with 
E >> 1. Then the averaged Jacobi equation for the limit e^ — )■ in the ultra- 
relativistic regime correspond to the differential equations 

/ 1 dX^^ dX"^ 

dXUX^ , dX^dX^ dXUiK 

For the above longitudinal electric field, the equations of motion are 

^ d^e" dc^ dx'' ^ , ^, de^ dx2 



,, /dX^ dX2 dX^ dX^ ^ . „, 

. , .z . .po ^ / dX^ dX^ dX^ d^ . 
d^e^ de° dX'' ^ , „, de'' dX^ ^ , „, 

^;^ /dX^' dX^ ^ . „, d^O dX^ ^ . „. \ 



dXi dX^ dX^di''. 



+ (''r%., + e'a,''r%.,)(^^ + 2 



d^X^ ,.n^1 .) ^ r,^2 MdX^ dX^ dX^ d^^ . 

In an inertial coordinate system the inertial terms are zero. Therefore, the sys- 
tem of equations in the linear longitudinal dynamics in the ultra-relativistic regime 
simplifies to 

^2^2 di^dX^^,^^ di^dX^^,^^ ,^AX^dX^^,^^ 

If e^ = {y^) -^ ss and since the distribution function has support on the unit 
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hyperboloid, {y'')-^ ^ 1 + a. Using also the decoupling condition -^^--^ ~ 0, the 
first two equations are 



Let us concentrate in equation (|6.2p . In the ultra-relativistic limit the velocity field 

~dr ~ 'y(*) (™ units where the speed of light is equal to 1). Therefore, the equation 
above can be written as 

"^^ +^E2-2jm'diE2{X + e) = 0. (6.3) 



di2 dt 



Let us approximate by Taylor's expansion E2{X+^^) = E2{X)+^^-^E2 in equation 



a 



()6.3p . Due to the translational invariance of the partial derivatives, di = ^ in the 
above expressions, by the chain rule. Then we have 

J2t2 j/:2 a 

^ + ^^2 - 2j{t)e^iE2{X + e) - E2{X)) = 0. 

If E2{X + S^) can be approximated linearly on ^, using Euler's theorem of homoge- 
neous functions one gets the following expression: 

^ + ^E2it) - 2^{t){E2{{X + - E2{X)) = 0. (6.4) 

This is the Jacobi version of the linear longitudinal dynamics. 

6.1 Examples 

1. Constant longitudinal electric field. 

In this case the equation of motion is 

(pe de ^ 

A particular solution is 

)^2 
^2 ^ _i_(-g-S2(t-to) _ 1)_ 

E2 

This equation is a linearized version of the longitudinal dynamics [lOj Vol.1, 
chapter 8]. 

2. Alternate longitudinal electric field. In this case, the electric field is of 
the form 

E2{X^ + f) = E2{0)sin{wrf{X^ + C^)). 
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The differential equation is 

-^+-j^E2iO)sin{wrfiX^+f))-2^E2iO){sin{wrf{X^H'^))-siniwrfX^)) = 0. 
We can expand this equation in ^, since ^ is smaU 

W2t2 Jt2 

-^+-^E2{0){sin{wrfX^)+cos{wrfX^)C^)-2-f{t) E2{0){cosiwrfX^)f) = 0. 
At first order in .^^ we have the equivalent expression 

(fi(:2 Jt2 

-^ + -j^E2{0)sin{wrfX^) - 27(t) E2mcos{wrfX^)e) = 0. 

We choose the initial phase such that sin{wrfX'^) ~ 0; therefore cos{wrfX^) ~ 
1 and the equation is 



£C 



2 



7 Discussion 



2j{t)E2m =0. (6.5) 



We have shown that the linear transversal beam dynamics in accelerator physics is 
obtained from the Jacobi equation of the averaged Lorentz connection for electro- 
magnetic fields F(X + ^) linear in the deviation variables ^. In particular we have 
proved that in the case when the magnetic fields are linear on ^, like in a dipole 
and quadrupole magnetic fields, the transverse dynamics can be interpreted as the 
dynamics of the Jacobi equation of an affine, symmetric connection ( ^V). A similar 
conclusion follows for the linear longitudinal beam dynamics, where the effect of 
constant and oscillating electric fields on the bunch of particles are described as an 
approximation of the Jacobi equation for the averaged Lorentz connection ('^V). 

One advantage of the theory presented in this work respect to the standard treat- 
ment of beam dynamics based on the one particle Lorentz force equation is that our 
theory involves only notions that are observable. In particular, the notion of refer- 
ence trajectory that we provide is linked with the mean velocity field of a kinetic 
model, which is an observable quantity. We also describe a relation between the 
distribution function and the collective behavior of the system. 

The function {Offc' ){t) measures the departure of the motion associated with 
the mean velocity vector field from the ideal reference trajectory associated with the 
trajectory of a point charged particle. 

Further theoretical research can be directed to generalize the approach to non- 
linear dynamics. This is apparently done by the so called generalized Jacobi equa- 
tions for affine connections [9]. Another direction is the study of conjugate points 
associated with the averaged Lorentz connection. Conjugate points can be asso- 
ciated with self-intersections of the solutions of the equation of motion. Finally, 
numerical implementation of the averaged Lorentz equation for numerical simula- 
tions can be interesting for accelerator particle physics, since its structure is more 
simple than the original Lorentz equation and it is expected better computational 
perspectives. 
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